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Abstract 



We prove weak-strong uniqueness results for the isentropic compressible Navier- 
Stokes system on the torus. In other words, we give conditions on a strong solution so 
that it is unique in a class of weak solutions. Known weak-strong uniqueness results 
are improved. Classical uniqueness results for this equation follow naturally. 



< 
rH ; 1 Introduction 



1.1 Presentation of the equation 



We shall study in this paper the Cauchy problem for the isentropic compressible Navier- 
Stokes system set in the d-dimensional torus T , with d > 2. The unknowns p and u stand 
^ I for the density and velocity of the fluid; they are respectively M"^, and M'^ valued and they 

Q,. ' are functions of the space variable x and of the time variable t. The system reads 

CS| ■ ( dtp + div{pu) = 

I->J: ! (1) < pdtu + pu ■ Vu - Lu + VP{p) 

O ; [ iP,u)\t=o = {po,uo) ■ 

^"^ ' where the notation u • V corresponds to the operator X^j^^ u^di. The pressure law P and 



^- (2) 



the Lame operator L are given by 

P{p) = p^ with 7 > 1 



L = ^A -h (A -F p)V div with /i > and \ + 2p>Q 



The pressure depends only on the density due to the isentropy assumption; the Lame 
operator arises since the fluid is assumed to be Newtonian. 

The mass 

■^{p) — pdx is constant. 

Formally, energy is also conserved 



£{p,u) =^ hlvWHml + -MtW, + r l|Vn(.)||i ds = hlVP^nog + -IIpoI 
2 1 Jo 2 7 



For the weak solutions presented in the next subsection however, it is not known whether 
the above energy equality holds, or simply the weaker energy inequality, that is the second 
equality sign above is replaced by <. 



The following transformation leaves the set of solutions invariant: 

(uo,Po) ^ {Xuo{Xx),po{Xx)) 

iu,p,P) ^ iXu{X'^t,Xx),p{X'^t,Xx),X'^P) . 

This is not a true scaling transformation since the pressure law P has to be scaled too. 
However if one disregards this problem and only focuses on u, one finds the classical scaling 
of the Navier-Stokes equation. Inspired by this, we will consider that a functional space for 
u and p to be scaling invariant if its norm is left unchanged by the above transformation. 

1.2 Weak and strong solutions 

The first line of research concerning solutions of ([T]) is strong solutions. Though this is no 
rigorous definition, we mean essentially by this solutions whose uniqueness (in the class 
where they are built up) can be proved ; they are generally constructed using iteration 
techniques. 

Danchin [3j built up global strong solutions for (small) data at the scaling of the equation. 
He had to assume that the initial density was bounded by below by a positive number. 
Cho, Choe and Kim were able to deal with data with vanishing density. However, 
they required a compatibility condition on their data and could not consider data at the 
regularity level of the scaling. 

Strong solutions have two major shortcomings: their local existence is not known below 
a certain regularity for the data, or in the presence of vacuum (vanishing density); and 
their global existence is known only under a smallness assumption. 

An explicit example of blow up has even been built up by Weigand [18j; but one has to 
apply an unbounded force to the system. 

Let us now examine weak solutions. These solutions are in general globally defined, built 
using a compactness argument, but their uniqueness is not known. 

Weak solutions for ([T]) have been proved to exist first under a spherical symmetry assump- 
tion by Hoff [9]. 

The case of general, finite energy data was treated by Lions [14j, with a technical restriction 
on 7, which was reduced by Feireisl [6] to 7 > 2- The weak solution which is obtained has 
finite energy, but not much more is known of its regularity. 

Jiang and Zhang [12j were able to build up weak solutions for any 7 > 1, for axisymmetric 
data. 

Finally, Desjardins [5] considered data whose regularity is slightly better than what finite 
energy would impose, {po,uo) £ L°° x H^. He could build up local weak solutions for which 
this level of regularity is conserved. However, their uniqueness is not known. Notice that 
L°° bound on po reminds one of the scale-invariant spaces discussed above. 
Of course, the basic problem with weak solutions is that the question of their uniqueness 
is not answered. 

1.3 Weak-strong uniqueness 

The idea of weak-strong uniqueness is the following: assume that a weak solution has 
some appropriate extra regularity (ie, is also a strong solution in some sense), and prove 
its uniqueness in the class of finite energy solutions, or some similar weak solution space. 



Weak-strong uniqueness gives in particular conditions under which the equation is well 

behaved, and weak solutions are unique. 

The first weak-strong uniqueness results were obtained by Prodi [16j and Serrin [17J for 

the (incompressible) Navier-Stokes equation, for which the same question can be asked. 

See the articles by the author [7] [8] for recent results about weak-strong uniqueness for 

incompressible and inhomogeneous incompressible Navier-Stokes. 

Weak-strong uniqueness criteria for ([T]) were given in [5]. 

2 Results obtained 

We present below our two theorems on weak-strong uniqueness for ([T]). 
It will be a convenient short-hand to denote 

Theorem 2.1 Suppose 7 > 2; o^nd take initial data {uq, po) such that 

po G L'^ and y/poUo G L^ . 
A solution (u, p) is unique on [0, T] in the set of solutions (p, u) such that 

(4) £{p,u)< 00 and VpeLpLy^^^) 

provided it satisfies, for a constant c, p > c> 0, and 

• ifd = 2and-f>2,Vue L\^L°° , LH G L^LP with p>2. 



• 



if d = 2 and -y < 2, Vu £ L\^L°^ , Lu G L^LP , Lu G L^''' L? with p > 2, 

ifd>3andj>2,Vue L\.L°° , Lu G L\L'^. 
• ifd>3and'y<2,Vu£ L\L°° , Lu G L|L^ , Lu G L^^'' L"^ . 

Remark first that this theorem can be extended to the case where the domain is the whole 
space R"' if d > 3. The proof applies then almost verbatim. 

Next we would like to mention a remark that was pointed out to us by Nader Masmoudi, 
and which yields a statement close to the above theorem. 

The idea is that the condition V/> G L^'^'Lv^^ '^' does not actually play any role in the 
estimates leading to the theorem, other than justifying formal manipulations. It therefore 
should be possible to try and dispense with it. A possibility for doing so is to include the 
inequality ([5]), which is the key point in the proof of the theorem, in the definition of a 
weak solution. To be more specific, we would define a new class of weak solutions, which 
would not only verify the energy inequality, but also inequality ([5]) for any p, u with the 
appropriate smoothness. It would be possible to adapt the argument used to build up 
global weak solutions in order to include this new condition, since inequality ([5]) holds for 
smooth p, u: this is the content of Lemma 13. 1[ 



So another version of the above theorem could be obtained, replacing the set given by (jl]) 
by the set of solutions (p, u) such that 

£{p,u) < cxD and {p,u) satisfies inequality ([5]) . 

A striking fact is that in this new version of Theorem 12. H weak strong uniqueness holds 
without any regularity requirement on the density. 

We now come to our second theorem. 

Theorem 2.2 Take initial data {uo,po) such that 

po £ L°° and ^/poUo S L^ ■ 

A solution (p, u) is unique on [0, T] in the set of solutions (p, u) such that 

p e LfL°^ ^pu G V^L^ Vu G L|l2 

provided that 

• ifd=2, 

Vp G L^L^ 

Vu G L^L°° 

Vi {dtu + u ■ Vu) G LtLP , 
with p > 2. 

• if d> 3, 

Vp G L^L'^ 

Vu G L^L°° 

Vt {dtu + u ■ Vu) G L^L'^ . 

Finally, it is interesting to compare the two theorems above... 

• ... but let us compare them first to existing results. 

Theorem 12.11 up to the technical condition Vp G LjJ L^^'< ''^ , only requires that 
the energy of {u, p) is finite, which is a very natural condition ; this is the first 
uniqueness result in such a wide class. Also notice that its proof is quite elementary. 
There are uniqueness results in classes of smoother functions : Danchin [3] and 
Hoff [lOj. A nice feature of the result of Hoff is that his criterion includes solutions 
that exhibit a particular kind of singularity, which is known to exist for (CINS): 
see Hoff [II]. 

As for Theorem 12. 2^ it is the first weak-strong uniqueness criterion at the scaling of 
the equation ; the result of Desjardins [5] is close, but misses the scaling. 

• Theorem 12.11 takes advantage of the whole energy (see the proof), whereas Theo- 
rem [22] mainly uses the finiteness of the kinetic energy ||y^u||2- For this reason, 
the first theorem above is not restricted to bounded densities. But the exponent 7 
is restricted to 7 > 2- Notice that this threshhold for 7 is also the one below which 
existence of weak solutions is not known. 



• In Theorem l2.2( one imposes that p is bounded. In some sense, it makes the pseudo- 
scahng transformation behave like a real scaling transformation (that is, it cancels 
the trouble with P) , and endows the equation with a real scaling. It is then no wonder 
that the conditions provided by Theorem 12. 21 correspond to this scaling (think of the 
case of Navier-Stokes, where the classical weak-strong uniqueness conditions are at 
the scaling of the equation) . 

• Finally, Theorem 12.21 includes the possibility of vanishing viscosity p. The estimates 
of Theorem 12.11 do not require that p stays away from zero, but making the ap- 
proximation argument work in this case would probably impose heavy additional 
technical restrictions, so we chose not to consider this case. 

3 Proof of Theorem [23] 

We consider two solutions of (CINS), {u, p) and (n, p), which share the same initial data, 
and are both assumed to be of finite energy. In the following, (n, p) will play the role of a 
general finite energy solution, whereas {u, p) will have additional regularity; our aim will 
be to prove the uniqueness of n, that is to prove that u = u. In order to achieve it, we 
shall estimate their difference with the help of the following lemma 

Lemma 3.1 Let {u,p) and {u,p) be two solutions of (CINS) verifying the assumptions 
of Theorem Ii2. 1\) . and 

def _ def 

U = u — u K = p — p . 



Set 

Fip,R) = -{R + pr-p'-'R--p 

7 7 

Then 

(5) 

dtWVpUWl + e||VC/||2 + -^dt\\F{p,R)h 

7-1 



;7 



< - f (p{U-Vu)-U + R^-U + jdivuF{p,R)j dx 



Remark 3.1 The quantity \\\/pU\\^-\ — ^||-F(/>, i?)||i is sometimes called relative entropy 
of p with respect to p. R seems to have been first used to prove weak-strong uniqueness 
results by Dafermos f^, who was studying conservation laws. More recently, Mellet and 
Vasseur f7^ . by making use of this quantity, obtained weak-strong uniqueness results for 
the one dimensional isentropic compressible Navier-Stokes equation. This relative entropy 
is useful in many other contexts : for instance, Berthelin and Vasseur JT^ relied on it to 
prove convergence of kinetic models to the isentropic compressible Euler equation. 

Proof of Lemma l3.lt We will prove the above Lemma without caring about regularity 
issues when manipulating expressions. A regularization procedure should be performed 
in order to prove it with full rigor. The assumption that u and u satisfy the hypotheses 
of Theorem 12.11 provides bounds on u and u that would enable one to carry out such a 
regularization step; we shall however skip this. 
First, subtracting the mass conservation equations for {u, p) and (n,p) gives 

(6) dtR + div(/>f/) + d\-v{Ru) = 



Subtracting the momentum conservation equations gives 

(7) {pdt + pu ■ V)U -LU + Vp^ - Vp^ = -R{dtu + u-Vu)-pU-Vu 

Substituting "~ ^ for dtu + u ■ Vu, this yields 



(8) 



{pdt + pu ■ V)U -LU + Vp"! - ^Vp^ = -pU ■ Vu - R^ . 

P P 



The next step is to take the (space) L^ scalar product of the above with U. Using the 
mass conservation equation, the first term gives 

(9) J{pdt + pu-V)U-Udx = ^dt\\VpUf2. 

The second term contributes, due to ^ 



(10) 



LU ■ U dx > e\\VU\\i 



with e > 0. The scalar product of U with the third and fourth summands of ([SD is a bit 
more tricky: 



(11) 



P. 
P 



Vp^ -^Vp' ] -Udx 



1 

7 



7- 



\.Udx = ^- I pUV (p^-i - p^'-^) dx 

— / div(pC/) {p'^~ — p^^ ) dx integrating by parts 

- I \dtR + div(i?n))(p^-i - f-^) dx by ^ 



7 



7-1 

7 
7-1 

7 
7-1 

7 
7-1' 



d 



dR 



a 



dtR-:—Fdx+ / u-VR-—Fdx+ / diYuR-—Fdx 



OR 



d 



dR 



dt Fdx - / dtp—F dx+ u- VF dx - u- "^PjpF dx + div uR—F dx 
dt I Fdx+ I div u (-F + p—F + R—Fj dx 

dt \ F dx + ^ I div uF dx . 



'27r(2^+i 



Notice that the assumption that V p as well as Vp belong to LjJ L^^i "^J implies p and 
p belong to L^ L'^"' . These two facts, combined with Vn , Vu G L'^L'^ ensure that all the 
expressions written above converge properly. 
Now putting together (f8l)-(f9l)- (fT0]) - (fTT]) . we get the desired estimate. ■ 



Lemma 3.2 With F defined as in Lemma \3.1[ for any 7 > 1 there exists a constant C 
such that 

R^<CF{p,R)p^-^ ifRe[-p,p] 
R"'<CF{p,R) ifR>p. 



Proof: It only consists of elementary computations, so we skip it. ■ 

Remark 3.2 Since p and p are non-negative functions, R = p — p is always larger than 
-P- 

As a last ingredient of the proof of Theorem I2.H we need the following lemma, almost 
identical to the one proved in [13j, page 78: 

Lemma 3.3 Take 7 > 2; '^^^ suppose that 



M.{p) = pdx >m and / p^ dx < M 

for positive constants m and M. Then there exists a constant C = C{m, M) such that 

\\v\\^^ < C iW^vh + WVvh) . 

Proof: Let us fix m and M and argue by contradiction. If the lemma is not true, there 
exists a sequence {pn^Vn) such that 

p'^dx>m>0 , pldx <M , \\^/p''v'^\\2 + ||Vt;„||2 ^ and ||v"|| -^ = 1 • 

We can extract a subsequence such that, for some a G M and p €z L'^, 

p^ ^ p weakly in U' and / Vndx ^ a . 

This implies 

/ 2d 
pdx >m and v" ^ a strongly m. L<i-^ . 

So a cannot be zero. Since 7 > 2 we deduce that 

^n|^n|2 _^ |^|2^ weakly in L^ , 

contradicting the fact that ||y^t'"||2 ^ 0. ■ 

Proof of Theorem 12.11 We present the proof only in the case where d > 3. The result 

2d 

is different in dimension 2 for then the Sobolev embedding H "^^ L^-a does not hold 
anymore; so one has to modify the proof in an obvious way. 

We will prove that the right hand side of ([5]) can be bounded by 

|l|V[/||i + /(t)(||V^t/||i + ||F||i) , 

with / G L^([0, T]); Theorem (j2.ip follows then by Gronwall's lemma. 

Bounding the first and third summands of the right-hand side of ([S]) is fairly easy: 



/ 



{p{U-Vu)-U + jdwuF{p,R)) dx < C7||V^||loo {W^Ug + \\F\\i) . 

7 



For the second summand, let us deal first with the region where R < p. One has then, by 
Lemma [321 R < C/F^^-t/^, hence 






Udx 



<cIvf 


Lu 

p-r/^ 


\U\dx 






< CIIFIll'^lll/ll ,, 

d — 2 


Lu 


by Hold( 

d 


<c\\F\\y^ 


i\\vu 


2 


+ \\V7^u\\2) 


Lu 

p7/2 



by Lemma 



<C 



Lu 



97/2 



IFWi + ^muh + W^uhY 



Let us now consider the region where R> p. 

If 7 > 2, by Lemma [3.2t the inequality R < CyFp^^'^''^ still holds, so we find, proceeding 

as above. 



R-^-Udx 
{R>p} P 



<C 



Lu 



97/2 



\F\\i + -^i\wu\\2 + \\^puhy 



If 1 < 7 < 2, we notice that R < p which gives, combined with Lemma 13.21 that R < 
/Fpi-'^/^ on {R> p}. Therefore 



Lu 
R^ ■ U dx 
{R>p} P 



< / s/F\Uy'-^p^-^'l'^\U\'^-^' 
<C\\F\\\'^\\Urjl\\^pUg-^ 



Lu 



dx 



Lu 



d 

7-1 



by Holder's inequality 



<C||F 



1 1/2 



y^uh + W^pUhV-'W^pUWt 



7-1 



,2-7 



Lu 



d 
7-1 



by Lemma 



<-^{\\yuh + \\^pUhY + c 



Lu 



3-7 



d 

7-1 



WWi + W^pUWl) . 



4 Proof of Theorem [22] 

We shall only consider the case d > 3; the case d = 2 can be proved following the same 

path. 

We consider a solution of ([I]), {p,u), satisfying the assumptions of the Theorem 12. 2 ^ and 

would like to prove its uniqueness in the class 



2 r2 



p G LfL°° y/^u G LfL^ Vu G L^L 

So let us pick (/J, u) a solution of ([1]) in the above class, sharing the same data as {p,u). 
In order to prove that (p, u) = {p, u) , we shall estimate the difference 

def _ def 

U = u — u R = p — p . 



Let us begin by estimating R. 



Lemma 4.1 The following inequality holds 

dt\\Rh < C ((IIpIIoo + ||/o||oo)||VC/||2 + \\Vu\U\R\\2 + ||Vp|U||?7||^ 

Proof: Subtracting the mass conservation equations for {p,u) and {p,u), we get 

dtR + pdivU + U -VR + u-VR + Rdivu + U ■Vp = . 

Taking the (space) scalar product with R and integrating by parts when needed gives the 
desired result. ■ 

Let us now estimate U. 

Lemma 4.2 If p and p are hounded by M in L°° , there exists C = C{M) such that for 
some e > 

i5t|IVP^Ili + e||VC/||i< ||i?||2||9t1i + n-Vu||d||C/||^+C'||V[/||2||i?||2 + ||^/pf/||il|Vn||oo • 

Proof: This inequality follows in a straightforward fashion after taking the (space) scalar 
product of d?]) with U and noticing that under the assumptions of the lemma, 

|(Vp^ - VpT , C/)| < C7||i?||2||Vf/||2 . ■ 

We also need to control the mean value of U. This is achieved by the following 
Lemma 4.3 The following inequality holds 

Udx < -— - (||p||oo||VC/||2 + ||Vn||2||ii||2) 

Proof: The lemma follows immediately from the following formula (appearing in Des- 
jardins [5]) 

and Holder and Sobolev inequalities. ■ 

As is always the case for this kind of problem, the proof of the theorem will be concluded 
by a differential inequality. The classical Gronwall lemma will not suffice here, so we will 
resort to the more general lemma proved in a previous paper by the author [8]. 

Lemma 4.4 Suppose that the following inequality holds 

(12) f'+{gf<af + (3g'g. 

where f , g' , a, and (3 are positive functions of the real variable such that 

f eL°^ , g{0) = , g' eL^ , a e L^ and Vi(3{t) G L^ . 

Then there holds 

■^«" -I -I fsP{s)'ds) f {gf + e-Soy(t) < /(O) . 







We can now proceed with the proof of the theorem. 

Proof of Theorem I2.2t We shall from now on assume that (p, u) and {p, u) satisfy the 

assumptions of the theorem. 

First, Lemma 14.31 and Sobolev inequality yield the control of \\U\ 



2d . 
d-2 



(13) ||f/||^ <C(||VC/||2 + ||Vn||2||i2||2) . 

d-2 

Inserting this in Lemma 14.11 gives 

dtWRh < C iWVUh + (||Vti||oo + \\Vuh) \\Rh) 
Integrating this last inequality, we obtain 



iRhit) <C [ \\VU\\2ds . 
Jo 



If we use this last inequality to estimate 11 -Rib and ()13p to control ||C/||^d_, Lemma 14.21 

d-2 

becomes 
1 



^dtWVpuWi + eWvuwi 

< C i \\dtu + u ■ Vn||d||VC/||2 / II VC/II2 ds + \\dtu + u ■ V-uHdH Vn||2 ( / || VC/||2 ds 

+ IIVC/II2 I ||V[/||2 ds + ||Vn||oo||^/p?7||i j 

<ci(\\dtu + u-Vu\\d + l)\\VU\\2 / ||VC/||2(is + t||atu + u- Vn||rf||Vu||2 ( / ||Vf/||^dsJ 

+ l|V^||oo|IV^?7||i) . 
In order to conclude, it suffices to apply Lemma 14.41 with 

fit) = \\\VpU[t)\\l + \ j'^WWUf^ds 
.'(*) = |livc/|b 

2 

a = 2 Vn 00 H — * Vu 2 9*^ + u ■ Vu\\d 
e 

P = -{\\dtu + u-Vu\U + l) .■ 
e 
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